The XENONIOO exclusion limit without considering iLgff as a nuisance parameter 
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In 2011, the XENONIOO experiment has set unprecedented constraints on dark matter-nucleon interactions, 
excluding dark matter candidates with masses down to 6 GeV if the corresponding cross section is larger than 
10^ cm . The dependence of the exclusion limit in terms of the scintillation efficiency (Les) has been debated 
at length. To overcome possible criticisms XENONIOO performed an analysis in which Xgff was considered as 
a nuisance parameter and its uncertainties were profiled out by using a Gaussian likelihood in which the mean 
value corresponds to the best fit X^ff value (smoothly extrapolated to zero below 3 keVnr). Although such a 
method seems fairly robust, it does not account for more extreme types of extrapolation nor does it enable to 
anticipate on how much the exclusion limit would vary if new data were to support a flat behaviour for Xgff 
below 3 keVnr, for example. Yet, such a question is crucial for light dark matter models which are close to the 
published XENONIOO limit. To answer this issue, we use a maximum Likelihood ratio analysis, as done by the 
XENONIOO collaboration, but do not consider Xeff as a nuisance parameter. Instead, Lgff is obtained directly 
from the fits to the data. This enables us to define frequentist confidence intervals by marginalising over Lf-g. 



I. INTRODUCTION 

After several decades of intensive searches, direct detec- 
tion experiments have now reached the required sensitivity 
to probe efficiently the parameter space associated with mas- 
sive weakly interacting massive particles (WIMPs). Just a few 
years after the world's best exclusion limits set by the EDEL- 
WEISS and CDMS experiments on the dark matter-nucleon 
elastic scattering cross section, the XENONIOO experiment 
demonstrated that the use of Xenon based technologies could 
actually beat Germanium detectors. By pushing down the ex- 
clusion limit by almost a factor 10 on the whole dark matter 
(DM) mass range (and with the present level of Krypton pu- 
rity 1 1 1), the XENONIOO experiment not only could rule out 
values of the DM-nucleon cross section as low as 7 10^^^ cm^ 
for DM particle masses of ^ 50 GeV at 90% confidence level 
Pi but could also exclude (similarly to CDMS) some of the 
light DM candidates which have been hypothesised to explain 
CoGeNT 12J, DAMA/LIBRA 13] , CRESST [4| and COUPP 
Q findings (unless one relaxes some assumptions as done in 
e.g. 1^-9 1 even though this may not be sufficient, e.g. 1 10 1). 

One element of controversy in the interpretation of these 
results is the dependence of this limit on the scintillation ef- 
ficiency of the detector (see for example ifTTU T?!). A differ- 
ent Xeff energy dependence at low nuclear recoil energy could 
indeed change the recoil energy associated with low mass 
WIMPs and possibly lead to a different exclusion limit than 
published in 1 1 1. To address this issue, the XENONIOO col- 
laboration used a profile likelihood analysis in which L^ff was 
taken to be a nuisance parameter and its uncertainties were 
profiled out with a Gaussian Likelihood 
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where tabs — represents the mean value of L^n smoothly ex- 
trapolated to zero (£eff = 0) at low energy and abs{t — tabs) = 
1, the la confidence region. A flat behaviour of L^ff at low 



energy would then be the upper limit of the 1-sigma contour.In 
principle, armed with such a modelling, the collaboration ac- 
counts for uncertainties in the extrapolation of Lgg at low en- 
ergy. This is important since the L^ff energy behaviour below 
3 keVnr is being currently debated. 

However, as the XENONIOO analysis is very complex and 
relies on many sources of uncertainties, the use of this method 
makes it difficult to determine what would be the real impact 
of a very different energy behaviour of Lgg (with respect to the 
mean considered by the XENONIOO collaboration) on the ex- 
clusion limit if the latter was determined from new data with 
unprecedented precision below 3 keVnr (or if some existing 
data were found to be less reliable than previously thought). 
This is particularly important in the context of light dark mat- 
ter candidates (eg. |14| ) where small recoil energies are ex- 
pected. 

Our work is essentially motivated by the fact that some light 
dark matter scenarios lie very close to the XENONIOO limit 
and should be therefore very sensitive to new measurements 
of Zeff below 3 keVnr. Several of these scenarios cannot be 
constrained by the LHC nor the measurement of the y-ray, 
synchrotron fluxes nor even directional detectors (see for ex- 
ample lil5nl8J ). Nuclear recoil direct detection experiments 
might be the only option to exclude or discover such scenar- 
ios, thus emphasising the need for a better determination of 
X-eff at low energy. In order to properly account for the lack of 
determination of the low energy behaviour of X-eff, it is neces- 
sary to make transparent the correspondence between the Lgff 
energy behaviour and the exclusion limit. Therefore, here, we 
will not treat Xeg as a nuisance parameter. Instead we will use 
directly the value of L^ff obtained from spline fits to the data 
and will not profile out the uncertainties (due in particular to 
the extrapolation at low energy) so as to quantify the effect 
Lgff has on the exclusion limit. 

The estimate of the uncertainties presented in this work is 
limited since we are not part of the collaboration. It is likely 
that we do not use up-to-date methods and data (as a matter 



of fact new data should be published soon). Also we base 
our analysis on several assumptions which are explained in 
the sections below. However, despite all these limitations, we 
could recover the exclusion limit that XENONIOO collabora- 
tion has obtained and can therefore highlight the large impact 
the low energy behaviour of Lf.f{ has on the exclusion limit. 
Note that in this paper we focus on Xgff uncertainties only; the 
astrophysical uncertainties will be addressed elsewhere. 

In Section III] we recall the spline interpolation to L^ff 
dataset as well as the extrapolation at low energies and discuss 
the robustness of the fit using an extended filter formalism. We 
use different types of interpolation and extrapolation (consis- 
tent with the 1-sigma contour defined by XENONIOO in [IJ). 
In Section [ni] we derive the exclusion limit for the mean Lgg 
interpolation and compute the exclusion limits for more ex- 
treme Xeff behaviour at low energies. Results are given in IV 
and conclusion inlVl 



II. Xeff 

The XENONIOO experiment aims at detecting dark matter 
particles via their elastic scattering interactions with Xenon 
nuclei in a two-phase (liquid and gas) time-projection cham- 
ber (TPC) detector A DM signal is then expected to have two 
signatures. The first one, referred to as the primary scintilla- 
tion signal ^i, arises directly from the interaction of a DM par- 
ticle with the liquid Xenon and measure the scintillation light 
in the liquid detector. The second, referred to as S2, happens 
in the upper part of the detector - at the liquid-gas interface - 
and measures the scintillation light which results from the drift 
of the free electrons that originate from the ionisation of the 
Xenon nuclei in the liquid phase after the DM interaction and 
which survived the recombination with ionised atoms. Both 
signals are measured in photon-electrons units (PE) |19| and 
are used to calibrate the detector's response to nuclear recoil 
events and ultimately to determine whether the experiment has 
actually detected dark matter events. 

The discrimination parameter is defined as 
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Events below the threshold of ER — —0.4 in the expected 
energy range are considered as potential DM events. The 
XENONIOO experiment uses the ratio of the two signals ^i 
and ^2 to discriminate between a DM and a background event 
so the identification of signal is actually sensitive to the pri- 
mary scintillation yield of recoiling Xenon nuclei in the liquid 
part of the detector As the measurement of the absolute scin- 
tillation yield is difficult, the quantity that is used by the col- 
laboration is the scintillation yield of nuclear recoils relative 
to that of 122 keV yrays from a ^''Co source. This is called 
the relative scintillation efficiency and is referred to as Lgg. 

The nuclear-recoil energy threshold E„r (in units of keVnr) 
of a signal is then determined by both Si and Zeff according 
to the relation, 
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where Ly is a normalisation factor for the light-yield of the 
122 keV gamma rays and Se, S„ are scintillation quenching 
factors for electronic and nuclear recoil respectively, due to 
the presence of an electric field (for XENONIOO, the values 
used are Se — 0.58 and S„ = 0.95). 

The determination of Ly and Zgff are therefore of utmost 
importance. While Ly has been measured precisely to Ly = 
2.20±0.09jjj^, there is no theoretical prediction for the en- 
ergy dependence of L^ff. An empirical formula was obtained 
in |20, 21] by fitting the data obtained in the same reference, 
namely 

Left — qnd qesc qel 

with q„ci the Lindhard factor (cf [201), q^sc reduction of the 
scintillation light yield and q^i a quench factor due to bi- 
excitonic collisions ll22l . 

Such an empirical fit reproduces the observation that the 
Lf.ff data decrease with decreasing energy and is also the as- 
sumption made by the XENONIOO collaboration in |[T] in or- 
der to obtain a conservative exclusion limit. Yet there are no 
measurement of L^s at low recoil energy. Besides, theoretical 
considerations by 1231 seem to favour a constant behaviour 
of L^ff at low energy. This would be consistent with the fit 
obtained by the XENON 10 collaboration |24| . 

The XENONIOO collaboration's strategy to incorporate the 
uncertainties on L^ff is to consider iLgff as a nuisance param- 
eter and profile out the uncertainties with a Gaussian Likeli- 
hood centred on the mean value of Lgg, that is the best fit. 
Similar assumptions are made for the other parameters which 
enter the analysis. Although this seems a robust approach, it 
is not very transparent. In particular, one loses the correspon- 
dence between the exclusion limit and the uncertainties on L^ff 
which arise due to a specific spline interpolation of the data 
and extrapolation at low energies. Indeed, the 1-sigma con- 
tour for ileff does not show on the exclusion curve obtained in 
im but since these uncertainties are due to the lack of data, one 
does expect to be able to keep track of them. In addition, it 
is hard to tell whether the final exclusion curve does take into 
account possible changes in the knots of the interpolation. 

In the following, we therefore adopt a different strategy. We 
still use a profile Likelihood analysis but we do not treat L^ff 
as a nuisance parameter. As a result, we can directly see the 
effect of the uncertainties on L^ff interpolation and extrapola- 
tion on the exclusion limit. We thus obtain several exclusion 
limits where the mean should be seen as the exclusion limit 
corresponding to the best fit of L^fi and where the edges of 
the contours correspond to the upper and lower parts of the 
//eff 1-sigma bands. Said differently, instead of obtaining one 
exclusion curve which would correspond to the best fit given 
all the uncertainties in the analysis, we prefer to draw the ex- 
clusion curves corresponding to the mean value and 1-sigma 
bands of £eff and let the reader marginalise 'by eyes' the ef- 
fect of Zeff on the exclusion curve. This approach enables us 
to anticipate the effect of a possible change in the physics of 
//eff below 3 keVnr 



A. X-eff interpolation 

To overcome the lack of knowledge about the low energy 
behaviour of iLgff, it was suggested by the XENON 100 collab- 
oration to perform an interpolation of the Chepel et al. li25J . 
Manzur et al. |20|, Plante et al. [26 1 and Aprile et al. 127] Leg 
data sets and perform an extrapolation below 3 keVnr Since 
older data sets (e.g. JSSl, ll29l,||30l,||3li) were disregarded in 
|[T], we will only consider them to understand their impact on 
the Zeff interpolation' . 

Like in LU, we perform a cubic spline interpolation to the 
four datasets previously mentioned and use five knots, placed 
at recoil energies of 5, 10,25,50 and 100 keVnr respectively. 
The best-fit cubic spline is found by freely varying the y-axis 
positions of these knots, while minimising the least-squares 
X^ goodness-of-fit parameter between the interpolated spline 
and the data (see lf33]| for a good discussion of the methodol- 
ogy). The result is shown in figurefl] along with the one sigma 
contour, obtained by looking for the maximum and minimum 
y-axis positions of the knots which satisfy x^ < x^,, + 5.89. 



positions (5 — > 10 keVnr and the additional knot at 75 keVnr) 
appears to be the enlargement of the errors in the extrapolated 
region for energies below the first knot. However there are 
also clear alterations to the interpolation around 75keVnr. 
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FIG. 1 : A fit of a natural cubic spline to data for the relative scin- 
tillation efficiency of Xenon, shown as a yellow line, along with the 
one sigma contour, shown in blue. The fit uses five knots, shown as 
red squares, at fixed positions on the x-axis of 5, 10,25,50 and 100 
keVnr. The uncertainty on the extrapolation is reflected in the top and 
bottom curves of the one sigma blue band. Note that recoil energy 
refers specifically to nuclear-recoils here. 

The choice of the x-positions of these five knots being 
somewhat arbitrary, we now perform another cubic spline in- 
terpolation where we place the knots at 10, 25, 50, 75 and 
100 keVnr. The translation of the lowest knot, from 5keVnr 
to lOkeVnr, has been performed to illustrate the effect of ig- 
noring the potentially less-reliable data below lOkeVnr. As 
can be seen in Fig. [2] the greatest change due to the new knot 



' An attempt was made by 1321 to measure I^if using the nuclear-recoil band 
of XENONIO. This data is not considered in our fits, but does provide 
an interesting alternative method of determining the relative scintillation 
efficiency of Xenon. 



FIG. 2: A fit of a natural cubic spline to data for the relative scin- 
tillation efficiency of Xenon, shown as a yellow line, along with the 
one sigma contour, shown in blue. The knots used to draw the best- 
fit spline are shown as red squares, at positions on the x-axis of 10, 
25, 50, 75 and 100 keVnr. The uncertainty on the extrapolation is 
reflected in the top and bottom curves of the one sigma blue band. 

Changing the knots influences the L^ff energy dependence. 
In particular, it changes the shape at high and very low energy. 
By adding a knot at 75 keVnr, we actually gave some weight 
to the single point at (55.2,0.268) which has for effect to drag 
the curve up around 50 keVnr. Removing the knot at 5 keVnr 
and instead extrapolating also changes the behaviour of L^ff 
at low energy. In particular, the uncertainties on Lgff become 
larger below 10 keVnr and notably the constant extrapolation 
moves to higher values of L^ff. 



B. Xeff extrapolation 

Since there are no data-points below nuclear-recoil energies 
of 3keVnr there is a great uncertainty on the energy depen- 
dence of L^ff at low recoil energies. The empirical behaviour 
found in EOl seems to imply that Lf.ff falls down to at low 
energy in a way which would be consistent with the spline fit 
of Zeff at higher energy. However, |23| suggests that based 
on the physics of Xenon recoil and an understanding of both 
the ionisation yield and scintillation efficiency, iLeff should be 
constant below 10 keVnr. Such an energy behaviour would be 
supported by |f33]| where it is argued that the drop in the scin- 
tillation efficiency observed by f20l could be due to the drop 
in sensitivity in the experiment. 

Given the lack of data, we will perform an extrapolation 
of our curves at low energy as in 1 1]. I.e. we adopt either a 
constant Zeff below a certain energy threshold or a drop to 0. 
For this latter case, we either extend the spline fit to 1 keVnr 
or to 2 keVnr (as in yj). The uncertainty on the extrapolation 



is reflected in the top and bottom curves of the one sigma blue 
band in both Figs{T]and|2] Finally, we also try a sharp cut-off 
of £eff at low energy for the bottom curve of Fig|2]in order to 
obtain the most conservative limit. 



C. Robustness of the fit 

Figures [T]andl2]show that even with slight modifications in 
the fitting procedure, the results for Lgg as a function of recoil 
energy can change significantly. In order to check the quality 
of a certain fit to the data, we employ the extended critical fil- 
ter formalism presented in 1 34 1 . This formalism finds a fit to a 
noisy data set by making use of the error statistics of the data 
points as well as a Gaussian prior probability distribution for 
the underlying curve. It is taking into account the possibility 
of outliers in the data, i.e. data points with significantly under- 
estimated error bars. This seems to be beneficial in the case 
of the Xeff measurements due to the wide spread and apparent 
inconsistency of the different data sets. 

Here, we feed the algorithm with different ileff-curves as 
mean for the Gaussian prior If the prior mean is akeady a 
sufficiently good fit to the data set, the result of the extended 
critical filter procedure will not deviate from it. If, on the 
other hand, the result of the data filtering differs from the prior 
mean input, it is a sign that the data prefer a different curve, 
even though the possibility of individual data points being out- 
liers is accounted for These outliers are accounted for in the 
algorithm by the inclusion of a correction factor for the er- 
ror bar of each data point (see f34l for all technical details). 
By narrowing the prior probability distribution for these cor- 
rection factors, we can force the algorithm to take each data 
point more seriously and thus find out which of the fits is most 
consistent with the data. 

In this way, we study the quality of the two cubic spline 
fits shown in Figs. [T] and [2] as well as the iLgff-curves given 
by the upper and lower one-sigma contours (i.e. the edges of 
the blue-shaded regions in Figs. [T] and |2]i. Using a reasonably 
wide prior for the error bar correction factors, we find that 
all of these curves are consistent with the data, except the top 
edge of the one-sigma region in Fig.fT] The exclusion of this 
one curve might, however, well be due to its behavior at large 
recoil energies and is likely not to be related to the extrapola- 
tion at lowest energies since the top one-sigma curve in Fig. [2] 
is not excluded although it is a more extreme extrapolation. 
Note also that the behavior at recoil energies below 3 keVnr 
is not constrained by this analysis. 

When narrowing the prior for the error bar correction fac- 
tors to more and more extreme shapes, more curves are suc- 
cessively excluded. It can thus be determined that the central 
fit in Fig. [T|is the most consistent one with the data. The mul- 
titude of Xeff-curves that is consistent with the present data, 
however, clearly underlines the importance of studying their 
influence on the resulting exclusion curve. In fact, yet an- 
other fit can be obtained by using a constant curve as prior 
mean for the extended critical filter and narrowing the prior 
for the error bar correction factors until deviations from this 
constant become significant. The resulting curve is shown in 



Fig. l3] along with the one-sigma contours of the two spline- 
fits shown in Figs. [T] and |2] 
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FIG. 3: Reconstruction of the calibration curve when using the ex- 
tended critical filter with a constant prior mean (corresponding to the 
mean of all data points), shown as a yellow line, along with the one 
sigma contours around the fits of Figs.[T|and[2] 



in. EXCLUSION LIMIT 

Now that we have determined the uncertainties on iLgff, we 
can compute the counting rate of dark matter events expected 
in the XENON 100 detector and deduce an exclusion limit for 
a given Xgff . For this purpose, we use a profile Likelihood ratio 
method and compute p-values for the signal and background, 
as done in |I] after randomly simulating 10000 'mock' data 
sets based on the XENON 100 data pubhshed in 1 1 J. 



A. Counting rate 

The recoil rate (per nucleus) is parameterised in the stan- 
dard form of (351, 
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where a is the WIMP-nucleus cross-section, q — \/2mi^E 
is the nuclear recoil momentum (with ota? being the nucleus 
mass), m is the WIMP mass, /j is the WIMP-nucleus reduced 
mass, p is the local WIMP density and r|(£',f) is the WIMP 
mean speed, given by the expression 
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In the above integral, u^it) is the relative velocity between the 
Earth-based detector and the WIMPs, with time-dependence 
arising from the motion of the Earth around the Sun, and 
Vn,in(£') is the minimum velocity for a WIMP producing a 
nuclear-recoil of energy E. Any astrophysical uncertainties. 



such as dependance on the galactic WIMP velocity distribu- 
tion /(v,Me(f)), will arise through this r|(£',r), and could affect 
the analysis. Note that the standard halo model is used in the 
proceeding analysis and we took the average velocity over the 
year and multiplied the final rate by 100 days (in accordance 
with H] where the total number of expected events assumes a 
run of 100 days with a fiducial volume of 48 kg). 

The WIMP-nucleus cross-section o can be further param- 
eterised (assuming spin-independent interactions, and equal 
coupling to protons and neutrons) as. 



a{E) = On 
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where On is the zero-momentum-transfer cross section 
for WIMP-nucleon interactions (hereafter any references to 
cross-section will be to On), A is the atomic mass, fjp is the 
WIMP-proton reduced mass and F{E) is the nuclear form 
factor (taken to be of the standard Helm form here) 11361 [37l . 
Equation ([TJ can now be used to calculate the signal rate per 
number of photoelectrons n in the detector: 
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where P{n,v{E)) refers to the Poisson distribution 
P{n,v{E)) = ^^ 138|. Here v{E) is the expected 
number of photoelectrons in the detector at an energy E and 
is given by the expression 



v{E)^ExUff{E)x^xLy, 
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where Lsf[{E) x Ly gives the number of expected scintilla- 
tion photons from a nuclear recoil per keVnr of nuclear-recoil 
energy E. 

In order to translate the expected rate per photoelectron ^ 
into something comparable with data, one must factor in the 
finite photomultiplier resolution OpM = 0.5 PE and knowledge 
about cut acceptance rjcuts- This gives an expression for the 
rate per nuclear-recoil signal (SI): 
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here G{Sl,n, \/napM) is a Gaussian distribution 
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It is then possible to calculate the expected number of signal 
events in the detector Ns, by integrating 



B. Profile Likelihood ratio 

The analysis described below follows the approach pre- 
sented in 1381 and centres around the application of the profile 
likelihood method. The analysis takes as input, both theoreti- 
cal parameters, such as the expected number of signal events 
Ns for a given WIMP mass m,^ and cross-section o, and a set 
of data-points. 



1. statistic test 

The so-called Profile Likelihood ratio for a particular 
dataset can be expressed as 
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where here L,„ax{o) refers to the likelihood function max- 
imised with respect to all parameters but a, which is held 
fixed, and Lmax{^) refers to the likelihood maximised with 
respect to all variable parameters^ respectively. The quantity 
q^, defined as 
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with a being the value of the cross-section which extremises 
the likelihood function measures the quality of the fit for that 
particular dataset and dark matter parameters (mass and cross 
section). The larger q^ is, the less signal-like these parameters 
are supposed to be, thus ruling out this particular value of a 
for a given dark matter mass and dataset. In principle, for 
every mass one should test the entire range of cross section 
that one initially considered. Here, however, we only tested 
values of the cross section which are relatively close to the 
XENONIOO Hmit by using a step-size of logio(o) = 0.02 in 
a range defined as [Oxemm, 10 O xenon] with Oxemm the value of 
the dark matter-nuclei cross section at the exclusion limit. 

The question of the set of data-points that should be con- 
sidered for the analysis is essential. In the following, we will 
use both the experimental dataset given in |JJ and simulated 
datasets that we will generate using a Monte Carlo. By consid- 
ering hypothetical alternate XENONIOO experiments, repre- 



sented by randomly simulated datasets (see subsectior III B 3 1, 
one therefore takes into account the random nature of the ex- 
perimental data. With this in mind, there should be a certain 
proportion of simulated datasets (the exact value of which de- 
pends on the chosen confidence) which provide a better fit 
than the actual experimental data. 

The log of the statistical test (qa) for a given dataset, dark 
matter mass and cross section is uniquely determined. We 
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these are specifically cross-section a, galactic escape velocity I'esc, total 
number of background events Nj; and the probabilities to be signal and 
background events Ej and E^ 



will refer to it as q'^^ when using the experimental data set 
and qa otherwise. To a given simulated data set corresponds a 
certain value of q^^ (for fixed dark matter parameters). Thus, 
one expects a q^ distribution, which can be used to define a 
p-value. 

The signal and background p-values {ps and py respec- 
tively), are defined as: 
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Here f{qa,Ha) is the probability density function (pdf) of all 
the qa values from simulated datasets, under the so-called sig- 
nal hypothesis //„, while f{qu,H()) is the pdf for the back- 



ground hypothesis Hq (cf subsection III B 3 i for fixed dark 
matter parameters. Here Hq refers to the electronic recoil 
events only while Ha refers to the electronic and nuclear recoil 
events. 

Under a desired confidence of 90%, one defines the exclu- 
sion curve by satisfying the condition that 
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for each dark matter mass. Incorporating the p-value for the 
background hypothesis Hq into the calculation of the exclu- 
sion curve enables to take into account the possibility that the 
background can mimic a WIMP discovery signal, by explic- 
itly requiring the chosen cross-section and dark matter mass to 
fit better to the signal hypothesis H^ than Hq. This is particu- 
larly important given the overabundance of background points 
compared to signal points for the XENON 100 data. 



The first term parameterises uncertainty in the escape ve- 
locity Vesc, which is treated as a nuisance parameter here. The 
second term compares, using a Poisson-distribution function, 
the number of expected data-points in each band /' for both 
signal and background, to the actual number of points iij. 
Here e^ is the probability for a background event to be in band 
j and Ei is the equivalent for signal events, as determined from 
calibration data. We start by defining e^ as e^ = inj^/Mi, where 
m^ is the number of electronic recoil data points in band j 
as appearing in the calibration data. Mi, = Jlj'^^i^ ^nd then 
marginalise over it. Similarly, before marginalisation, el = '^ 
where mi is the number of nuclear recoil data points in band j 
as appearing in the calibration data Ms = Y.j'^i- The product 
elNs represents the expected number of signal events which 
fall into the nuclear recoil bands (given a particular cross sec- 
tion, WIMP mass and choice of /-eff). 

Since the expected total number of signal events Ns is a 
function of cross-section a, data-points in bands where ei > 
e^ will have the greatest effect on the best-fit a value for a 
particular dataset. In particular the analysis is very sensitive 
to any data-points appearing in the lower bands (see figures [5] 
andpl, where electronic recoil/background events are unlikely 
to occur. 

Finally the last two terms parameterise the uncertainty in 
probabilities el and e^, due to the expected Pois son- variance 
of the number of ^'^'AmBe and ^''Co calibration data-points 
in each band, ml and m^ (here M, and M/, are the total num- 
bers of nuclear and electronic recoil points respectively for the 
calibration data). 



2. Likelihood Function 

Our likelihood function is the same as in I.38J except that 
we do not parameterise the uncertainty in the relative scintil- 
lation efficiency so as to make explicit the impact of £eff on 
the exclusion limit. 

Like ||38]| . we make use of bands to discriminate between 
electronic and nuclear recoils in S1-S2 space, thereby sepa- 
rating the signal from the background, allowing a more strin- 
gent limit to be placed on WIMP mass and interaction cross- 
section. 

The likelihood function (L) is given by the equation below: 

-i^ ^ ./v(,*^'obsj^escj X 
/ 23 



lYlPinJ,ieiNs + eiNb))x 
P{mlieiMt))xP{mi,ielMs)) 



where here /' is run over each of the 23 bands (the bands them- 
selves are shown in figure l5]l, and rtj is the number of data 
points in band j. 



3. Dataset Simulation 

Since simulated datasets play a vital role in determining the 
pdfs f{qa,Ha) and f{qo,Ho), it is important to discuss their 
method of generation, and the uncertainties involved. 

As discussed in the previous section, by defining an exclu- 
sion curve at 90% (or any value different from 100%) , the nat- 
urally random nature of the experiment is taken into account. 
Since the XENON 100 experiment can only ever be performed 
once, it is possible that any observation, or non-observation, of 
possible signal data-points could be due, wholly or in-part, to 
statistical fluctuations. One seeks to improve this possibility 
by positing hypothetical alternate XENON 100 experiments, 
which differ only in their sampling of the statistics. Practi- 
cally these alternate data-sets are represented by simulations, 
based on the actual experimental data [IJ. 

Since the simulated data should attempt to mimic that ob- 
tained by the experiment, the data-points must be arranged on 
the S1-S2 plane, where SI is equal to the number of photo- 
electrons observed for a particular event, and S2 represents 
the ionisation yield. This can be achieved with knowledge of 
the expected distribution of nuclear-recoils, associated with 
WIMP signal events, and electronic-recoils, associated with 
background events, in S1-S2 space. Such information is con- 



tained in the calibration data obtained by the XENON 100 ex- 
periment, who used a ^''Co source for samples of electronic 
recoils and an ^^' AmBe source for nuclear recoils. 

In practice this calibration data is binned and normalised, 
to give a pdf for signal and background events, shown in fig- 
ure |4] Given a desired number of signal (nuclear-recoil) and 
background (electronic-recoil) events, a simple Monte Carlo 
algorithm can be used to generate simulated datasets. 




FIG. 4: Contour plots of calibration data on the S 1 -S2 plane, showing 
the distribution of nuclear recoils in red and electronic recoils in blue. 



It is possible to extend the simulation algorithm further, and 
improve the accuracy of the simulated datasets, by incorporat- 
ing information about the WIMP energy spectrum, equation 
(|6]l, into the determination of the nuclear-recoil data-points^. 
However, since there are very few candidate signal points seen 
in the data, and to avoid possible problems with bias, this ex- 
tension has not been incorporated into the current analysis. 
Even so, the uncertainty in exactly how to simulate datasets 
most accurately will contribute a source of uncertainty to the 
final exclusion curve. 

The above method was used to generate 10000 simulated 
datasets, with an expected number of 2 signal events (nuclear- 
recoils) and 534 background events (electronic-recoils), be- 
tween 51=4 and SI = 20, as seen in the data obtained by 
the run of the XENON 100 experiment 1 1 1 after 100 live days 
of data-taking ^. In addition to these 'signal + background' 
datasets, so-called 'background-only' datasets were gener- 
ated, with an expected number of 536 background events and 
no signal events. 



^ Such an extension is not necessary for tlie background points, as electronic 
recoils have been shown to have a flat spectrum at the energies considered 
here 1 19] 

'' As an aside, by extending the SI region down to zero, discrimination be- 
tween signal and background points would be easier in principle, due to 
the clear difference in shape between the probability contours of nuclear 
and electronic recoils at low SI, as seen in figurel4] However the increase 
in sensitivity is granted at the cost of greater susceptibility to systematics, 
especially the il^ff uncertainty. 



A plot of one such signal-i-background simulated dataset is 
shown in figure |5] The analysis itself is blind to whether a 
point was generated as a nuclear or electronic recoil, how- 
ever the fitting of the cross-section is more sensitive to the 
lower bands, where fewer background events are expected. 
Due to the abundance of electronic-recoil events compared to 
nuclear-recoils, determining which points are due to which is 
a difficult challenge, and so a clearer discrimination between 
signal and background only arises statistically when consider- 
ing many such datasets, motivating the choice of a confidence 
limit other than 100% for the XENONIOO limit. Hence, even 
with high statistics, the ability of the analytical tools to dis- 
criminate signal from background is limited, contributing a 
natural source of error to any determination of the best-fit val- 
ues of the cross-section and number of background events, 
and so ultimately to the final exclusion curve. 




FIG. 5: An example of a simulated dataset, with two nuclear-recoil 
(signal) events, shown in red. The rest of the points are electronic- 
recoil (background), shown in blue. The black lines divide the S1-S2 
plane into the bands used for the analysis. 

Values of q^ were calculated for each dataset under the 
prescription of section |IIIB| The signal-i-background and 
background-only q^, values were then binned separately into 
two normalised histograms (for each value of cross-section 
and WIMP mass), to give the pdfs /((7o,//o) and f{qa,HQ) 
respectively. In this way, the signal-i-background datasets rep- 
resent the signal hypothesis H^,, as they are generated under 
the assumption that the two candidates-signal events seen in 
the XENONIOO data are in fact due to nuclear recoils. Con- 
versely the background-only datasets take these points to be 
due to background electronic-recoils, thereby coming under 
the background hypothesis. 

An example of f{qa,Ha), for a specific WIMP mass and 
cross-section, is shown in figure |6] Note that, due to Wilks' 
theorem one expects the pdf to approach a %^ distribution as 
the number of sampled datasets increases, a trend which is 
indeed observed. Each pdf is fitted with an analytical x^ func- 
tion, to speed up computation and to avoid susceptibility to 
statistical fluctuations at higher values of qu. The value of a 
which satisfies equation ( [T3| l is sensitive to this fit, for both 





WIMP Mass (GeV) 



FIG. 6: The probability distribution function of q^ under the sig- 
nal hypothesis, for a WIMP mass of m = 7 GeV and o^f = 2.51 x 
10^ cm". The pdf has been constructed using 10000 simulated 
datasets, assuming 2 expected signal events and 534 expected back- 
ground events. The actual points of the pdf are shown as bars, while 
the best-fit y} distribution is shown as a black line. The blue dashed 
line indicates the value of ^^ for the real experimental dataset. 



FIG. 7: The XENONIOO exclusion curve using the best-fit cubic- 
spline for Xeff from figure [T] shown in red, along with the one sigma 
systematic uncertainty due to Xgff from the fit of figurefT] in blue, and 
the one sigma uncertainty from the fit of figure[2] in yellow. Note that 
all exclusion curves have a natural uncertainty of logjg o = ±0.02. 



f{qa,Ha) and f{qa,Ho), and so any uncertainty in the best- 
fit x^ function will contribute to the uncertainty in the final 
exclusion curve. 



IV. RESULTS 
A. Systematic Error In XENONIOO Exclusion Curves 

The profile likelihood analysis has been performed using 
the best cubic spline fit to iLgff, along with the top and bot- 
tom edges of the one-sigma bands from figuresfl] and |2] The 
resulting exclusion curves are shown in figure |7J 

Clearly the systematic uncertainty due to the relative scin- 
tillation efficiency is appreciably large for WIMP masses be- 
low 10 GeV, with the majority of the variation arising from 
the extrapolation-uncertainty for L^ff at low nuclear-recoil en- 
ergies. The lower edges of the uncertainty bounds on figurelT] 
coiTespond to the upper edges of the one-sigma regions on fig- 
ures[T]andl2] hence a flat extrapolation of Zeff at low-energies 
tends to result in a stronger XENONIOO exclusion limit. 

The upper curve in Fig |8] represents the most conservative 
exclusion limit that one can derive from the present L^ff data. 
This illustrates that the conclusion from the XENONIOO col- 
laboration is fairly robust above ^^20 GeV and uncertainties 
at large mass are really small, as claimed by |[T1. 

In comparison to similar works such as fT2l, which anal- 
ysed the first set of XENONIOO data without using the Profile 
Likelihood method, our conclusions are similar, but not iden- 
tical, in the low-mass region of parameter space, for the case 
where their cut-off of events at 51 = 1 was relaxed. However, 
in the current analysis it is possible that the use of the Profile 
Likelihood method has changed the size of the systematic L^ff 




WIMP Mass (GeV) 



FIG. 8: The XENONIOO exclusion curve as shown above with also 
the expectations if Lgff has a sharp cut-off below 10 keV The red 
plain curve is the mean fit; the blue dotted curve represents the 1- 
sigma contour corresponding to the bottom curve in Fig. [2] and the 
green dashed line represents the exclusion curve associated with a 
Lgff function with a sharp cut-off below 10 keV. 



errors for low-masses due to increased sensitivity, in addition 
to the different choices of spline interpolation for L^-ff itself, 
compared with lfT2l . 

It should be noted that there are a variety of uncertainties 
affecting the exclusion curves from the analysis, primaiily the 
uncertainty in fitting the x^ distribution to the pdfs of q^, but 
also the flexibility in the actual method of dataset simulation, 
and the overlap of background and signal regions on the S1-S2 



plane (as discussed in section III B 3 1 
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B. Discussion and Implications 

The relative size of the variation of the XENON 100 exclu- 
sion curve with L^s at low masses compared to that at high 
masses, can be understood in terms of the WIMP recoil spec- 
trum, an example of which is shown in figurel9] and equations 
(|4|to(|6]). 



larger n. Conversely for larger masses, all values of n con- 
tribute terms of the same order to (|6]), hence there will be no 
dramatic change in the final result, due to the relative scintil- 
lation efficiency. 

Additionally, in the case of the forms of Zgff from figures 
[T]andl2] the systematic uncertainty of the exclusion curve for 
low WIMP masses is further amplified, relative to the higher 
masses, due to the larger uncertainties of L^ff at low energies. 




1.0 1.5 2.0 

Recoil Energy (keV) 



FIG. 9: The recoil spectrum for a WIMP of mass 5 GeV and cross- 
section of Oat = lO^-'^cm^. From zero to 2.5keVnr the differential 
reaction rate gj changes by many orders of magnitude. 

The equation for ^ Q has two terms in the integrand: the 
WIMP recoil spectrum ^ and a Poisson distribution. The 
Poisson term is peaked at a particular value of energy, which 
increases for larger numbers of photoelectrons n. Hence for 
a certain value of n, there will be a region along the energy 
axis, of g|, which contributes most to the integral. By chang- 
ing the functional form of L^s, the value of the nuclear recoil 
energy at which the Poisson distribution peaks will change 
(by approximately 1 keVnr for the different parameterisations 
considered here). Hence for a particular value of n, the inte- 
gral of equation (|4]) will now receive dominant contributions 
from different areas of ^ along the energy axis, when iLgff is 
altered. 

To explain the different variation in the final exclusion 
curve due to iLgff seen at large and small WIMP masses, one 
must compare recoil spectra. For low masses, the recoil spec- 
trum changes rapidly at low energies (see e.g. figure l9]), be- 
fore falling off at a few keVnr. However, at larger masses, the 
spectrum is largely constant until much higher energies be- 
fore falling off; in the case of m^^ — 500 GeV the cut-off is at 
approximately lOOkeVnr. Clearly, the greatest change of ^ 
with //gff will be seen for values of n where the Poisson term 
of equation Q is peaked at recoil energies where ^ varies 
most rapidly. Since the peak of the Poisson in energy is at 
larger values for higher n, the largest variation in ^ will be 
seen for higher masses at high n, while for lower masses it 
will be seen predominantly at low values of n. 

Finally, equation (J6| contains a sum of ^ over all n. For 
low WIMP masses, the lower values of «, where the great- 
est variation due to Lgff occurs, dominate over the terms with 



V. CONCLUSION 

In this paper we have assessed the uncertainties on the 
XENON 100 exclusion curve 1 1 1 due to the lack of knowledge 
about the low energy behaviour of the scintillation efficiency 
of Liquid Xenon detector Our analysis is motivated by the 
existence of low mass dark matter scenarios (below 10 GeV) 
which lie close to the published XENON 100 limit (see for ex- 
ample |IT-T7l). 

The use of a profile likeUhood analysis (in which uncertain- 
ties on Lgff were profiled out) enabled the XENON 100 collab- 
oration to obtain an exclusion limit which is free from large 
uncertainties [1 J. In particular, the limit on low mass WIMPs 
seems very precise while one would expect to recover at least 
the 1-sigma uncertainty band which accounts for the lack of 
determination of Xeff below 3 keVnr 

In order to understand this behaviour, we have performed 
a similar profile likelihood analysis but did not consider Lgff 
as a nuisance parameter Instead L^ff is defined directly from 
the fits to the data. We show that the exclusion limit obtained 
by the XENON 100 collaboration at high energy is very ro- 
bust. The uncertainties on the exclusion limit due to Xeff are 
very small and all our exclusion curves are very similar to the 
XENON 100 exclusion limit. Such a conclusion was not nec- 
essarily obvious since different types of interpolation of L^ff 
data give different behaviours for Left, even at high energy. 
However this can be understood from our robustness of the 
fit analysis which shows that all spline fits at high energy are 
equally good. Hence heavy dark matter scenarios close to the 
XENON 100 Umits will not be affected by a better determi- 
nation of Xeff at high recoil energies. This implies that dark 
matter scenarios just above the XENON 100 limit are probably 
excluded indeed (provided that other sources of uncertainties 
which are not accounted for in this analysis are not too large) 
and, those just below, very close to be ruled out. 

At low recoil energies^, our results (cf FiglTJi show that 
the mean value of Xeff (which is in agreement with the mean 
Lf.ff curve considered by the XENON 100 collaboration) gives 
an exclusion limit that is similar to (although stronger than) 
the XENON 100 limit but a more extreme behaviour of iLgff 
at low energy (cf the flat extrapolation or a sharp cut-off be- 
low 3 keVnr for example) leads to a very different exclusion 



^ Our analysis of the robustness of the fit shows that the data still prefer 
a smooth cut-off to zero for ileff below 3 keVnr, even though a different 
behaviour has been suggested from theoretical arguments I23l . 
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limit. Should new data favour such a type of behaviour for 
Leff at very low energy with unprecedented precision, the ex- 
clusion limit would be different from that presented by the 
XENON 100 collaboration, even though such a behaviour was 
included in the 1- and 2-sigma contours considered for Zeff by 
the XENON 100 collaboration. 

Perhaps a reason for this discrepancy is the Gaussian Like- 
lihood term which was assumed for Xeff in the XENON 100 
analysis (and which was centred on the mean value of L^ff) 
together with flat priors. Combined with the other Likelihood 
terms which describe the uncertainties from the analysis itself, 
it may be that the XENON 100 maximum likelihood analysis 
over-fits Zeff and thereby suppresses alternative possibilities 
for this crucial quantity, leading to over confidence in exclud- 
ing light DM scenarios. 

In any case, since various iLeff fits give different exclusion 
limits, it seems more conservative to "track" the effect of 
different Xeff energy behaviour at low energy on the exclu- 
sion curve and define (frequentists) confidence intervals by 
marginalising over the different fits to Zgff data. With this 
method, one is in principle ensured not to bias the analysis 
towards the present best fit curve of Lgg since it may in fact 
not be the correct function to consider given the lack of data 
in this energy region. Finally, we note that a flat behaviour 
of Lgff (as suggested by 1231 ) at low recoil energy would ac- 
tually set a stronger exclusion limit. Hence the need for new 
measurements of L^ff at low energy. 

Our analysis was based on a likelihood chosen so as to keep 



the current analysis as close as possible to the one done by the 
XENONIOO collaboration ||I][T9l|38|, however it should be 
possible to improve the method in identifying WIMP signals 
over background. If the S 1 -S2 plane were to be divided into a 
grid of bins, then, on top of the knowledge of the actual struc- 
ture of the set of data-points in this space from calibration 
data, one could use additional knowledge of the recoil spectra 
for background and signal. Such a spectrum is already known 
along the SI -axis for WlMPs (equation (|6])), and is distinct to 
that from electronic recoils, which is constant for low energies 
|fT9]| . This would provide an additional method of discrimina- 
tion, which would be especially important if a future run of 
the XENONIOO experiment were to claim a discovery signal, 
but is difficult without more knowledge of the experimental 
set-up. Additionally, it should be possible, in principle, to 
test the efficiency of a particular likelihood function in recon- 
structing theoretical parameters such as the cross-section from 
given simulated datasets, since the parameters used to gener- 
ate them are known (see ||39]| for a similar discussion). 
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